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Abstract 

The thermal fields of two Bridgman-like configurations, representative of real systems used 
in prior experiments for the detached growth of CdTe and Ge crystals, are studied. These de- 
tailed heat transfer computations are performed using the CrysMAS code and expand upon our 
previous analyses [14] that posited a new mechanism involving the thermal field and meniscus 
position to explain stable conditions for dewetted Bridgman growth. Computational results 
indicate that heat transfer conditions that led to successful detached growth in both of these 
systems are in accordance with our prior assertion, namely that the prevention of crystal reat- 
tachment to the crucible wall requires the avoidance of any undercooling of the melt meniscus 
during the growth run. Significantly, relatively simple process modifications that promote fa- 
vorable thermal conditions for detached growth may overcome detrimental factors associated 
with meniscus shape and crucible wetting. Thus, these ideas may be important to advance the 
practice of detached growth for many materials. 
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1 Introduction 


The detached Bridgman growth is a new technique of crystal growth without crucible contact. 
The dewetting phenomenon was observed for the first time in microgravity experiments [1,2] then 
successfully reproduced in earth experiments [3-10]. The prevailing wisdom for conducting detached 
growth on earth is to establish a liquid meniscus between the crystal and the crucible wall and 
maintain it by controlling the pressure difference between the gas space above the melt and the gap. 
This practice follows early theoretical analyses of the stability of the process that concentrated on 
meniscus mechanics and indicated that dewetted growth depends on the wetting angles of the melt 
with the ampoule and the crystal and the pressure applied along the gas-liquid meniscus [11-13]. 

However, many of the experiments conducted at favorable theoretical growth conditions accord- 
ing to meniscus mechanics have shown reattachment events during the growth process. In a prior 
thermal-capillary analysis of detached Bridgman growth [14], we have put forth a new explana- 
tion for process stability that invokes heat transfer arguments to explain reattachment events that 
would end dewetting. According to our prior results, crystal reattachment to the crucible wall is 
explained by the freezing of the liquid meniscus between the growing crystal and the crucible wall, 
which occurs when the solidification isotherm passes above the liquid-gas surface of the meniscus. 
This work also identified possible remedies to prevent freezing of the meniscus, including changing 
the meniscus shape via wetting angle changes and increasing the pressure under the meniscus (al- 
though this approach is limited by the possibility of meniscus rupture by gas bubbling). A more 
robust approach is the tailoring of thermal conditions to coax the solidification isotherm to pass 
through the gap, under the meniscus, rather than through the melt, above the meniscus. Such 
conditions were shown in [14] to correspond with the promotion of convex solid-liquid interface 
shapes. 

The present work is meant to extend the understanding provided from our prior analysis [14] 
by focussing on thermal conditions associated with furnace design. We carry out detailed heat 
transfer computations for two configurations like those employed in prior experimental studies of 
the detached growth of germanium by Palosz et al. [9] and cadmium telluride by Chevalier et al. [5]. 
Our analyses reinforce the notion that favorable thermal conditions may prevent reattachment 
during detached growth. 

2 Model 

Our prior analysis employed a rigorous, thermal-capillary model that concentrated on the detailed 
coupling of interfacial and field phenomena but applied a very simplified representation of furnace 
heat transfer [14]. Here, we apply the sophisticated CrysMAS code developed by the Crystal 
Growth Laboratory of the Fraunhofer Institute of Integrated Systems and Device Technology (USB) 
in Erlangen, Germany [15] to perform rigorous simulation of heat transfer in multi-zone Bridgman 
furnaces. While this promises a much greater degree of realism with respect to modeling furnace 
heat transfer, the drawback of this approach is that we do not fully represent the coupling of 
meniscus dynamics with the system. Rather, here we simply prescribe the position and shape of 
the meniscus based on our prior results [14]; more details are provided below. Neither our prior 
modeling approach [14] nor the one taken here is entirely correct; however, both reveal interesting 
behavior and useful understanding of these complex systems. 

We apply CrysMAS for the global modeling of the heat transfer in the whole system, including 
the ampoule and contents, the furnace, and thermal shields. We consider heat transfer occurring 
by conduction and radiation in an axisymmetric configuration. Since our prior analysis indicated 
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that convective heat transport in the melt was unimportant in these small-scale systems [14], we 
choose to ignore melt convection in the current model. We approximate the effects of gas-phase 
convection in the furnace by employing heat transfer coefficients and effective thermal conductivities 
in open regions. Radiant heat transfer is computed through all gas-filled regions (including the 
gaps between the crystal and crucible and between crucible and ampoule) using view factors and 
enclosure analysis. Unless specifically stated otherwise, radiant transport is assumed to occur 
through all quartz components and is computed using a semitransparent band approximation, 
where the material is assumed to be transparent for radiation wavelengths of 1-4 /im and opaque 
for all other wavelengths. A thermal shield, assumed to have a thermal conductivity k — 20 
Wm -1 K -1 and an emissivity e = 0.8, is considered to surround each furnace in our simulations, 
and a constant ambient temperature of T — 300 K is imposed as an outer boundary condition. 
Additional furnace details will be discussed in the results sections. 

We employ the powerful inverse simulation mode of CrysMAS to compute input power to 
resistive heaters in order to satisfy temperature control set points. We also apply this approach to 
impose the constraint of the solidification temperature to the trijunction point. This allows for the 
joining of the prescribed fixed-radius crystal, the assumed meniscus shape, and the computed (and 
varying) solidification interface in a mathematically and physically consistent manner. 

Thus, computations proceed in the following manner. A small gap of constant thickness 0.1 mm 
is applied between the crystal and the inner crucible. Then the vertical location of the trijunction, 
where gas, melt, and crystal meet, is specified. As mentioned above, the meniscus shape is not 
computed by the code but is set by drawing a curve connecting the trijunction and the crucible 
wall, so that the wetting condition at the wall is attained. Finally, inverse simulation is applied to 
compute the power to the heaters in order to reach the solidification temperature at the trijunction, 
which defines the outer terminus of the solid-liquid interface and, in turn, sets the form of the 
thermal field of the entire system. 

3 Results 

3.1 Two- heater Bridgman furnace 

Figure 1 shows the simplified geometry used to model the two-heater Bridgman furnace employed 
for the detached growth of the germanium by Palosz et al. [9]. The geometry is considered to be 
axisymmetric, with the system centerline drawn on the left side of the schematic. The outer silica 
ampoule has an inner tube containing the charge. Several sets of experiments were conducted in 
this configuration using different inner tubes, ostensibly to investigate different wetting conditions 
on detached growth. The furnace has two heating elements separated by an adiabatic zone. To 
drive directional solidification, the furnace heaters are translated upward at a constant rate. 

Palosz et al. [9] found that successful dewetting over the entire growth run only occurred when 
a pyrolytic boron nitride (BN) inner tube was employed and attributed this success to the high 
value of the Ge contact angle on BN surfaces ( 6 ~ 170°). Other crucibles proved less successful in 
attaining stable, detached growth conditions. 

3.1.1 Effect of meniscus shape 

We first consider the conditions of successful detached growth of Ge in BN crucibles. As shown in 
Figure 2, the specific geometry consists of a BN inner crucible within a quartz ampoule; the gap 
between the crucible and ampoule is filled with gaseous argon. The position of the trijunction is 
set to be exactly in the middle of the adiabatic zone by adjusting heater power according to the 
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procedures described in Section 2. We concentrate here and in subsequent computations of this 
system only on the relative position of the solidification isotherm, which is critical to ascertain 
conditions that may avoid reattachment events [14]. 

In the detailed view shown on the lower left-hand side of Figure 2, we consider a meniscus shape 
representative of the real system, with a wetting angle of 9 ^ 170° and a convex- concave shape 
that is comparable to those computed by the thermal-capillary model in [14]. The thermal field in 
this case shows that the solidification isotherm is very closely aligned with the liquid-gas surface 
of the meniscus. Although there is a very small region of undercooled melt near the wall wetting 
point, most of the meniscus is kept at a temperature above the melting point, and this situation is 
favorable for dewetted growth, as observed in experiments. (Note: as will be seen in Section 3.1.3, 
the situation becomes more favorable as growth proceeds.) 

As a thought experiment, a second computation was carried out with a different meniscus shape 
but keeping all other parameters identical to the first computation. Here, a much decreased gap 
pressure and a smaller wetting angle are considered to produce a concave meniscus shape, as seen 
in the upper detailed view of Figure 2. In this case, the solidification isotherm is positioned above 
the meniscus, giving rise to a large undercooled zone in the melt. This is expected to quickly 
solidify, and the crystal will be attached to the crucible. These results demonstrate the importance 
of the meniscus geometry, as determined by pressure and contact angle, and are consistent with 
prior stability ideas based on meniscus mechanics [8,16-19]. However, the source of instability in 
this computation has its origins in heat transfer, as argued in [14]. 

3.1.2 Effect of ampoule-zone alignment 

Based on the ideas first put forth in [14], we undertook a series of computations to ascertain 
thermal conditions in this furnace that would enhance the stability of detached growth, especially 
for meniscus shapes that are less favorable than that obtained in the highly wettable BN crucible. 
Interesting results are achieved by simply changing heater settings to push the position of the solid- 
liquid interface upward. Figure 3 shows the predicted outcomes for an interface position raised 3.5 
cm above the adiabatic zone and into the hot zone of the furnace. 

Insets on the left-hand side of Figure 3 show the relative positions of the meniscus and solidifica- 
tion isotherm for three different crucible configurations. The upper case represents the BN crucible 
considered in the prior section. Under the new conditions, the solidification isotherm still tracks 
closely with the meniscus shape; however, the small undercooled melt region observed in the prior 
case (see Figure 2) is now gone. The middle case represents an inner quartz crucible held tightly 
within the outer quartz ampoule; the meniscus is assumed to form a straight line and wets the wall 
at a smaller angle of 6 ^ 160°. As is evident from the relative position of the solidification isotherm, 
this configuration presents a more favorable thermal condition than the BN case, with the isotherm 
clearly passing below the meniscus. A yet better situation is presented when the space between 
the inner crucible and the outer ampoule is filled with gas having a low thermal conductivity, as 
shown in the bottom inset. Here, even the worst-case scenario of a very low pressure with a concave 
meniscus shape still produces a thermally stable configuration. These results imply that favorable 
thermal conditions, here obtained by selected ampoule-zone alignment, may promote conditions 
to achieve and maintain detachment, even for mechanically stable meniscus shapes that are less 
conducive to detachment. 
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3.1.3 Stable growth 

Finally, we carry out several sets of time-dependent computations that move the heating zones 
upward to simulate conditions for growth. The inverse simulation for the initial condition supplies 
power levels for the heaters, which are subsequently held constant as the heating zones are moved 
and transient results are computed using forward simulation. 

Our first set of computations was performed for the BN crucible configuration described in 
Section 3.1.1, with an initial condition corresponding to the BN case shown in Figure 3. We do 
not show the outcome of this series of computations but state that stable thermal conditions were 
obtained for all stages of growth and that the estimated ratio of the solidification growth rate, V g , 
to the furnace translation rate, Vt, was V g /Vt ~ 1.5. This compares very favorably to the observed 
experimental results [9]. 

We next performed a series of thought experiments designed to test the assertion of the prior 
section that favorable thermal conditions may be able to compensate for less favorable wetting 
characteristics. Namely, we consider the outcome of growth using the ampoule-zone alignment of 
the prior section with quartz-crucible configurations employed by Palosz et al. [9] which were prone 
to reattachment events in their experiments. 

The first of these thought experiments considers the quartz crucible and ampoule with a gap 
that was discussed in the previous section. Figure 4(a) shows the initial condition for growth, with 
the interface positioned at 3.5 cm above the adiabatic zone, corresponding to a seed length of about 
2 cm. In this initial state, the solidification isotherm passes through the gaseous gap, as shown in 
the detailed view. The thermal conditions become more favorable for detached growth when the 
furnace is translated up, as is shown in Figure 4(b) and (c). In this system, the solidification rate 
was slightly greater than the BN case, with V g /Vt ~ 1 . 66 . 

The second thought experiment is depicted in Figure 5. Here, we add a graphite layer between 
the inner quartz crucible and the outer quartz ampoule. The major effect of the graphite layer 
replacing the air gap of the prior case is to move the position of the solidification isotherm upward 
due to the greater amount of heat transfered through the high-conductivity graphite. Though less 
favorable than the previous case, this system still exhibits thermal stability against reattachment 
during the entire growth run. Another effect arising from the graphite-promoted heat transport 
through the crucible-ampoule assembly is that the growth rate is quite similar to the furnace 
translation rate in this system, with V g /Vt ~ 1 . 1 . 

3.2 Three-heater vertical gradient freeze furnace 

We next consider an idealized vertical gradient freeze (VGF) configuration with three heaters, shown 
in Figure 6 and representative of that used for the detached growth of CdTe by Chevalier et al. [5]. 
As in the prior model, the geometry is considered to be axisymmetric, and the system centerline 
is drawn on the left side of the schematic. Directional solidification is obtained by powering down 
the furnace heaters to achieve a constant rate of temperature reduction at control set points. 

Chevalier et al. [5] attained successful detached growth only for an optimized thermal field in 
the furnace. The CdTe charge, with dimensions 1.4 x 8 cm, was contained in a silica crucible 
contained by an outer silica ampoule. The furnace had three heating elements put at a distance 
of 1 mm from the outer ampoule. The temperature inside the furnace was measured by using 
thermocouples placed in the space between the crucible and ampoule. According to Chevalier 
et al. [5], the optimal set of temperatures measured in the experiment at the beginning of the 
growth process was T\ = 1, 163 K, T 2 = 1, 188 K, and T 3 = 1, 373 K. This set of temperatures has 
been imposed at three control points (see Figure 6 ), and the power supplied in the furnaces has 
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been computed by using inverse simulation. 


3.2.1 Effect of meniscus shape 

We first compute the thermal fields at the beginning of the solidification process, shown in Fig- 
ure 7, for the optimal case reported in [5]. The position of the solid-liquid interface, calculated 
from inverse simulation requiring the thermal field to meet the optimal set-point temperatures given 
previously, corresponds to a seed length of approximately 3.7 cm, which is in good agreement with 
the experimental report of a seed of 3.5 cm [5]. The interface is concave with a small deflection 
of approximately 0.6 mm. The upper detail shows a case where the meniscus shape is set accord- 
ing to our prior thermal-capillary computations for the growth of CdTe in a quartz crucible [14] 
corresponding to a wetting angle of 9 = 100° and an applied pressure difference of A p = — 10 
mbar. Under these conditions, the solidification isotherm passes through the gap, just below the 
meniscus, therefore thermal conditions are favorable to support dewetting at this initial stage of 
growth. The lower case shows a computation under the same conditions, where the meniscus shape 
has been simplified as a straight line. Since there is virtually no change to the thermal conditions, 
the straight-line meniscus shape is used in ensuing computations for this system. 

3.2.2 Effect of crucible heat transfer 

Knowing that dewetted growth in this system was particularly prone to reattachment events for 
most furnace thermal settings, we consider another series of thought experiments where we examine 
how thermal stability conditions may be affected by crucible design. In Figure 8, we consider a 
later stage of growth, corresponding to 5.9 cm length of the solid, for several different crucible 
situations. 

Figure 8(a) shows the system considered in the prior section that featured a semitransparent 
quartz crucible. Even though the initial state shown in Figure 7 was stable, crystal reattachment 
becomes likely by the time the crystal has grown to the height shown in Figure 8(a), as indicated 
by the thin undercooled liquid layer along the entire meniscus. This somewhat puzzling result is 
at odds with the reported successful growth experiment of Chevalier et al. [5]. We therefore seek 
an explanation by examining the details of crucible heat transfer. 

To bracket the possible heat transfer behaviors produced by different crucibles, we first consider 
the case of a quartz ampoule that is totally transparent to radiation, shown as Figure 8(b). Here, 
stronger cooling caused by increased heat transfer through the crucible makes the undercooled melt 
layer significantly expand. This situation is less likely to lead to successful growth than the case of 
the semitransparent crucible. 

Better results are obtained when we consider a semitransparent quartz ampoule whose exterior is 
coated with a layer of opaque graphite, exhibited in Figure 8(c). In this system, thermal conditions 
move the solidification isotherm downward, so that it passes through the gap rather than the melt. 
Finally, Figure 8(d) presents a case where the quartz ampoule is considered to be opaque to thermal 
radiation. This case also offers a stable thermal configuration for maintaining detached growth, 
with the solidification isotherm passing beneath the melt meniscus. Both of these cases produce 
outcomes that are more consistent with the reported successful case of Chevalier et al. [5] . 

3.2.3 Stable growth 

The thermal field for the opaque crucible case of the prior section (Figure 8c) is very similar to that 
of the graphite-coated, semitransparent crucible (Figure 8d). Since the opaque-crucibe case is also 
far less computationally taxing, we choose to employ it in transient computations for a growth run. 
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We compute the course of a growth experiment by reducing the temperature of the control points in 
time. Figure 9(a) show the initial condition, and Figures 9(b) and (c) show times that correspond 
to a lowering of temperature of 60 K and 80 K, respectively. Clearly seen in the detailed views is 
that the solidification isotherm remains always below the meniscus so that the entire growth run 
exhibits stable conditions against reattachment. 

Interestingly, the behavior of this case is very similar to those of the transient Bridgman simula- 
tions (see Section 3.1.2), especially with respect to how the thermal field becomes more favorable as 
the solidification interface position moves higher within the hottest zone of the furnace. Accompa- 
nying this more favorable thermal field for dewetted growth is the evolution of a convex solid-liquid 
interface shape. 

4 Conclusions 

The global modeling of the thermal field in two vertical Bridgman-like crystal growth configura- 
tions has been performed to better understand favorable thermal conditions for successful detached 
growth of Ge and CdTe crystals. From a crystal growth modeling perspective, the model employed 
previously in [14] was a first step toward a fully integrated thermal-capillary analysis, as has been 
implemented for other meniscus-defined growth systems [20-26]. The analysis employed here rep- 
resents another step toward a model that more faithfully depicts the details of heat transfer in tech- 
nically realistic furnaces. We believe that the next advance will effectively couple these two models 
with their respective strengths using algorithms that are currently under development [27-30] . 

The computations for the specific systems considered here show some interesting results. We 
simulated the vertical Bridgman growth of Ge by Palosz et al. [9] and found, consistent with their 
results, that the large wetting angle of germanium on boron nitride surfaces was an important 
factor to promote conditions for stable, detached growth. Surprisingly, our computations predicted 
that simply reorienting the ampoule with respect to the furnace zones substantially improved the 
thermal conditions needed to avoid reattachment. By initiating growth much higher into their hot 
zone, favorable thermal conditions were predicted for continued detachment, even for systems that 
did not exhibit large wetting angles. In addition, substantial improvements in heat transfer, that 
affected the placement of the solidification isotherm with respect to the meniscus and hence the 
prospects for successful detached growth, were attained by relatively simple changes in crucible- 
ampoule design. 

Analysis of the vertical gradient freeze configuration developed by Chevalier et al. [5] for the 
detached growth of CdTe yielded initially puzzling results. Namely, our computations predicted 
that, despite an initial condition that would support dewetted growth, continued growth under 
reported conditions would eventually lead to undercooling of the meniscus. These poor stability 
characteristics were eliminated when the transparency of the inner crucible was altered, either by 
a graphite coating or by making it opaque to radiation. Hence, we speculate that the success 
reported in [5] may have arisen from intentionally applied graphitization of the crucible (although 
such procedures were not reported) or by some source of surface contamination of the crucible 
that deposited an opaque coating. Improved thermal conditions were also predicted, for coated or 
opaque crucibles, when the solid-liquid interface moved higher into the hot zone during growth. 
Thus, the conditions needed to maintain dewetted growth in this system improved in a manner 
similar to that in the Bridgman system when the growth interface was placed deeper into the hot 
zone. For both systems, thermal conditions that promoted a convex solid-liquid interface shape 
also produced favorable conditions to maintain detached growth. 
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Our prior work [14] argued that the long-standing notion of dewetted growth stability being 
determined solely by meniscus mechanics is incomplete and that thermal effects are also profoundly 
important. The results presented here strengthen that argument and propose important ideas 
for advancing the practice of detached growth. Clearly, the properties of the meniscus are very 
important in these systems. First, the meniscus must be mechanically stable to even allow the 
detached geometry to be established and maintained. Second, large wetting angles between the 
melt and crucible wall and large applied pressures in the gap are beneficial, since they move the 
melt adjacent to the meniscus upward and potentially away from the solidification isotherm [14]. 
What is new about the understanding obtained here is that the effects of less favorable meniscus 
shapes, such as those that are set by smaller wetting angles or lower applied pressures, that lead to 
reattachment under nominal conditions, might be ameliorated by the engineering of favorable heat 
transfer conditions in the system. Though supported only by computational evidence at this point, 
we propose that relatively simple modifications to current practice may attain thermal conditions 
allowing contactless growth during the entire process. Thus, these ideas may be important to 
advance the practice of detached growth for many materials. 
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Figure 1: Schematic diagram of the computational domains for the Bridgman configuration with 
two heaters employed by Palosz et al [9]. 
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Figure 2: Relative positions of the solidification isotherm and meniscus predicted for the configu- 
ration using a BN inner tube; the position of the solid-liquid interface is located in the center of 
the adiabatic zone. Two detailed views show cases for a convex-concave meniscus (below) and a 
concave meniscus (above). 
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Figure 3: Relative positions of the solidification isotherm and meniscus predicted for several dif- 
ferent ampoule-crucible configurations; the interface is positioned at 3.5 cm above the adiabatic 
zone. 
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Figure 4: Evolution of the position of the solidification isotherm when the furnace is translated up 
to simulate growth; the ampoule configuration consists of two quartz tubes separated by gas: (a) 
initial state; (b) heaters translated by 2 cm; (c) heaters translated by 3 cm. 
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Figure 5: Evolution of the position of the solidification isotherm when the furnace is translated up 
to simulate growth; the ampoule configuration consists of two quartz tubes with a layer of graphite 
between: (a) initial state; (b) heaters translated by 2 cm; (c) heaters translated by 3 cm. 
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Figure 6: Schematic diagram of the computational domains for the VGF configuration with three 
heaters employed by Chevalier et al. [5]. 
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Figure 7: Thermal field predicted for the initial configuration of successful detached growth de- 
scribed in [5]. Two detailed views show relative positions of the solidification isotherm and meniscus 
for a convex meniscus (above) and a straight-line approximation to the meniscus (below). 
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Figure 8: Thermal field predicted for an advanced position of the interface, corresponding to 
5.9 cm length of the solid . Detailed views show relative positions of the solidification isotherm 
and meniscus. Case for (a) semitransparent quartz ampoule; (b) completely transparent quartz 
ampoule; (c) semitransparent quartz ampoule with its exterior coated by a graphite layer; (d) 
opaque quartz ampoule. 
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c) 


Figure 9: Evolution of the position of the solidification isotherm when the furnace is powered down 
to simulate growth under the conditions described in [5]: (a) initial state (seed length of 3.7 cm); 
(b) solid length of 5.9 cm; (c) solid length of 7.1 cm. 
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